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ABSTRACT 

Context. Photon orbital angular momentum (POAM) is normally invoked in a quantum mechanical context. It can, however, also be 
adapted to the classical regime, which includes observational astronomy. 

Aims. I explain why POAM quantities are excellent metrics for describing the end-to-end behavior of astronomical systems. To 
demonstrate their utility, I calculate POAM probabilities and torques from holography measurements of EVLA antenna surfaces. 
Methods. With previously defined concepts and calculi, I present generic expressions for POAM spectra, total POAM, torque spectra, 
and total torque in the image plane. I extend these functional forms to describe the specific POAM behavior of both single telescopes 
and interferometers. 

Results. POAM probabilities of spatially uncorrelated astronomical sources are symmetric in quantum number Such objects thus 
have zero intrinsic total POAM on the celestial sphere, which means that the total POAM in the image plane is identical to the total 
torque induced by aberrations within propagation media and instrumentation. The total torque can be divided into source- independent 
and dependent components, and the latter can be written in teiins of three illustrative forms. For interferometers, complications arise 
from discrete sampling of synthesized apertures, but they can be overcome. POAM also manifests itself in the apodization of each 
telescope in an array. Holography measurements of EVLA antennas observing a point source indicate that ~ 10% of photons in the n 
= state are torqued to n t states. 

Conclusions. POAM quantities represent excellent metrics for characterizing instruments because they are based on real physics and 
are used to simultaneously describe amplitude and phase aberrations. In contrast, Zemike polynomials are just solutions of a differ- 
ential equation that happen to ~ correspond to specific types of aberrations (e.g., tip-tilt, focus, etc.) and are typically employed to fit 
only phases. Possible future studies include forming POAM quantities with real interferometry visibility data, modeling instrumen- 
tal aberrations and turbulence of the troposphere/ionosphere in terms of POAM, POAM-based imaging algorithms and constraints, 
POAM-based super resolution imaging, and POAM observations of astrophysically important sources. 

Key words, instrumentation: interferometers - methods: analytical - methods: data analysis - techniques: image processing — 
techniques: interferometers — telescopes 

1. Introduction 



lEliasI (l2008l) developed extensive semi-classical formalisms to describe photon orbital angular momentum (POAM) in astronomy. 

He assumed spatially incoherent sources, so these formalisms are significantly different from those used in laboratory situations. He 

concentrated more on instrumentation rather than astrophy sics, and included first principles, concepts, definitions, calculi, examples, 
l^ and applications. 

rS As a general rule, aberrations that damage wavefronts and POAM spectra should be minimized in both hardware and software: 

jrt \ Primum non torquer^ The amount of damage to wavefronts and POAM spectra should be estimated and removed by off-line 

■ ■ ■ image processing algorithms whenever possible. I now extend previous POAM work toward imaging metiics for single telescopes 

and interferometers, leading to a deeper understanding of how such instruments work. 

2. Basic Concepts 

The apertures of all real astronomical telescopes are finite in size, so they are only capable of producing diffraction-limited images 
up to a resolution of ~ D"' {D is the aperture diameter expressed in wavelengths). Assuming that the aperture response is uniform 
or at least azimuthally symmetric (radially apodized), this loss of information manifests itself as blurring. The only rigorous way to 
reduce blurring is to increase the aperture size. 



** The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 

' This Latin phrase - loosely translated as "Above all, apply no torque" - is a shameless adaptation of the most widely quoted words from the 
Hippocratic Oath "Primum non nocere," which means "Above all, do no harm." 
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Conversely, if the aperture response is azimuthally asymmetric (or, azimuthally apodized), the image is subject to non-uniform 
distortion. These distortions may occur over all angular scales larger than ~ D"'. The only rigorous way to reduce distortion is to 
eliminate the asymmetry in the aperture response. 

I distinguish between these two types of image degradation because it is much harder to build a very large telescope compared to 
figuring a very smooth mirror of nominal size. There is yet another reason for making this distinction, namely that radial apodization 
does not modulate POAM spectra, while azimuthal apodizatio n does modulate POAM spectra. In other words, aberrations modulate 
each input POAM state into one or more output POAM states (Elias 2008) and introduce image distortion. 

Interferometers measure visibilities at discrete points in the synthesized aperture. In turn, these measurements are transformed 
into "dirty" images of astronomical sources. This process is mathematically equivalent to punching small pinholes into the opaque 
aperture of a large single telescope and imaging the resulting interference optomechanically. 

The act of aperture sampling is a form of azimuthal apodization that modifies POAM spectra. Other physical effects - including 
instrumental imperfections, turbulence in the troposphere or ionosphere, etc. - manifest themselves as sample amplitude and phase 
eiTors and also modify POAM quantities. Symmetrically reducing the weighting of long baselines, on the other hand, only blurs an 
image and doesn't modulate POAM quantities. 

I can press the analogies even further. Image-processing algorithms, such as CLEAN and MEM, remove artifacts to yield a 
model of the true source up to a certain resolution. These artifacts are produced by asymmetric apodization, which means that 
image-processing algorithms actually estimate and eliminate changes in POAM spectra. 

3. Definitions 

According to lEIiasI (l2008h . the total POAM of a wavefront arising from the celestial sphere, in units of fi, is 

Iz = Lzlfi = 2_i "^Pf^.m^ (la) 

m=-oo 

where 

OS 

Pm,m — ^m,m I ^ ~ ^m,m I ? ^m,m \^^) 

m=-oo 

is the probability that a single photon (or the fraction of many photons) is in POAM state m, Bm,m is the (m, m)''' POAM autocorre- 
lation, and B is the total intensity. The ensemble of probabilities represents the intrinsic source POAM spectrum. The total intensity, 
which is integrated over the celestial sphere, is identical to the sum over all POAM autocorrelations. The same formulae can be used 
to describe the total POAM of a wavefront incident upon the image plane by replacing Iz — » fe, Lz — > Lz, Pm,m — > Pm.m^ B^^n — * 
B,„„„ and B ^> B. 

I define the total torque as the difference of the total POAM between the image plane and celestial sphere, or 

DO oa 

m=-oo m=-oo 

where the ensemble of t„, „ comprise the torque spectrum. Torque is normally defined as change in angular momentum per unit 
time. Since "per unit time" is ambiguous in this context I ignore it, so torque has the same units as POAM. 

The electric fields of a "natural light" astronomical source projected onto the celestial sphere are spatially uncorrelated. I call 
this scenario the "Standard Astronomical Assumption" (SAA), which led to the creation of the POAM calculi (Elias 2008; Tables 
1-4). Because of SAA, £-,„-„, - B„„, and p-m-m = Pm.m for all m, implying that lz = 0.1 derive this result in Appendix lAl Equation 
|2]then becomes 

oo oo 

m=— oo m=—o3 

Torque applied by propagation media and instrumentation modifies the source POAM spectral components as they travel to the 
image plane, but the total POAM in the image plane does not depend on the total POAM from the celestial sphere because of the 
symmetry in m. Although this equation is independent of lz, the + 1 transitional probabilities of the source do affect lz (cf. Section 

nil except for point sources at the center of the field of view (FOV). 

Maser photons traveling through turbulent gas and photons scattering off" Kerr black holes (iHarwi t 2003; T amburini et al. 11201 Ih 
may not satisfy the SAA condition. Therefore, these sources can exhibit Iz + Q and are beyond the scope of this paper. 

4. Single Telescopes 

I derive POAM quantities for single telescopes first bec ause it is rel atively simple to do and the results can be extended to other 
types of instruments (e.g., interferometers; cf. SectionlSTl. lEliasI (l2008l) created generic calculi that can describe the POAM response 
of any instrument, so I will use them here. 
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4.1. Initial Mathmematics 

Consider Figure [1] the schematic diagram of a single telescope looking at an object on the celestial sphere. The intensity response, 
in system form, for an SAA source is given by 



Bin) 



E(n ; t): 



I' 



; d^nD(Q.M)E(n;t) 

2/1^2 

2 



= fd^nD(ii,n) l^ Eiiiu) j= fd^Q.P(nM)Bin), 



(4a) 



where il = (pcos0,psin(/i) is the coordinate on the celestial sphere, n = (p'cos0',p' sin(^') is the coordinate in the image plane, 

B(il ) and B{ii) are the intensity distributions (E(il ; /) and £(£1; /) are the corresponding electric fields), P{il , fl) = f (Q - £1) is 
the point-spread function (PSF), 



Din , n) = D(n - n) 



'!"'' 



-j2!rRiil -ii) 



<R) 



(4b) 



is the diffraction functon (DP), R = (R cos ifr, R sin t//) is the coordinate in the aperture (in units of wavelength), and j'(R) is the func- 
tional description of the aperture apodization (cf. Section|2]i. Any optical system, including propagation media and instrumentation, 
that can be exp r essed in this mathematical form can be expanded into any of the POAM calculi. 

From lEliasI (l2008t Row 3 of Table 3), the (m, m)''' single-telescope SAA POAM autocoiTelation density is 



Bm.mip') = i-\E,„(p' 

where 



->l1 ^ / 



d^^p,n,,n(p',n)B(n), 



(5a) 



EJp'-t) ^ iz \ d(f>'e-J""^' E{n;t) ^ j d^Q.D„,(fi',n)E(n-t) o £(Q ; f) = ^ £„,(p'; 



1 
27r Jo 



t) 



Jmip' 



(5b) 



is the m image-plane POAM state. 



Pm,m(p', n) 



D,n{p',n) 



Yj f^m,k{p\p)e 



-jkij, 



k=-oo 



is the (m, m)''' PSF sensitivity. 



1 



2n 



D„(p',n) = -^ I d(f>'e--'""i'' D(n,n) o D(n,n) ^ V D,„(p',n)eJ""^' 

2^ Jo „,tl. 



(5c) 



(5d) 



is the m DF sensitivity, 



ni<lel 

Jp,gip' ,P) ^2jt I dRR JpilnRp') Sp-giR) Jq{2nRp) 
Jo 



is the (p, q) integral function, R,ci is the telescope radius, Jg(x) is the g order Bessel function of the first kind, and 

-2n 



1 r -> T -> 

Sg(R) = — dif,e-J^^'<'siR) o s(R) = V SgiR)e^«'^ 

2^ Jo „M, 



(5e) 



(5f) 



is the g'^ azimuthal Fourier component of the aperture function s(R). Integrating Equation |5a] over radius in the image plane leads 
to the {m, m)''^ POAM autocorrelation 



Bm,m = lim 2;r 

Prov— >°o 



XPfov n 

dp' p' Bm^,n{p') = I 



Bm,m{p') = #Qn,,m(n)B(n), 



where pfov is the FOV of the image plane, and 



Prav->°o 



2;r dp'p'P^^mip'M) ^ > > /-'27r dRRs^^tiR) sl_i{R) Jk{2nRp)J,(2nRp) 



-j(k-t)^ 



(6a) 



(6b) 



is the (m, m)''' PSF sensitivity kernel. I assume that pfov is large enough to capture most of the radiation scattered through the 
telescope aperture into the image plane. The complete derivation of Equations |5]|6] may be found in AppendixlBJ 
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All POAM quantities defined in Section[3]can be formed from Equations|6h-b. For example, the total POAM in the image plane 



IS 



1°° \ r \ °° ^l->ir ->^ 

fe = -g 2 '«Bm,m = -g M^Q Yj '"^'«."'(") Bin) = -g I d^n£,(n)B(n), 



(7) 



where Xz(fi) is the total POAM kernel. Total POAM on the celestial sphere for all SAA sources is identically zero, or /z = (cf. 
Appendix [All. This statement means that the total POAM measured in the image plane is identical to total torque applied to the 

wavefronts. Therefore, I interchangably employ the quantities Iz <-* t and -Cz(fi) <-* T'ifl), where T(il) is the total torque kernel. 
For a more physical understanding of Equation|7] Equations |6^-b are combined with various mathematical identities to create three 
"illustrative forms" of Iz which emphasize different aspects of POAM for single telescopes (cf. Section l4~2T l. 

4.2. Illustrative Forms oflz 
The first illustrative form of Iz is 



Iz^ Y^ mp1„,„ + lm27T I dppln \ dRR (ZnRp) 



Y Pn,n±l(p) Y P'km^l^'^^ 



(8) 



where pJJ, „, is the m''' POAM state probability in the aperture, Pn,n±i{p) is the transitional probability density between POAM states 
n and n + 1 on the celestial sphere, and pl^^^^^iiR) is the transitional probability density between POAM states m and m + 1 in the 
aperture. The transitional probabilities correspond to ±1 selection rules. I derive these equations and the variables contained therein 
in Appendix ICl 

The second illustrative form of Iz is 



'-Z 



mp„ 



pco pR 

I dpp 2n I 
Jo Jo 



Im27r I dppln 



dRR (InRp) 



iS^xip)' 
B 


S^iiR) 
S 



(9a) 



where the Sq:i(p) are the first order rancors of the source, the S±i(R) are the first order rancor sensitivit ies, and S i s the integrated 
squared magnitude of the aperture function. This expression proves that rancors and rancor sensitivities (lEliasll2008l) . 



S„(p) = — j d<t,e-JP^B{il)^ Y ^"'.'«-» « B{il) = 2 !Bp{p)eJP^ 

^■^ ^ 111 — — d-vi n — , d-v^ 



Bi,k(p) = {-Ei(p-t)El{p-t) 



(9b) 



Spi^R)^:^ I dipe-iP'^Si^)^ ^ S„^,,n-p{R) « 5(R) = ^ Sp{R)e'P^ 

*^^' m——oo p=—oo 



Si,k(.R) = Si(R)sl{R), S(R) 



^(R) 



(9c) 



calculated directly from intensities and squared aperture functions (instead of electric fields and aperture functions) and related to 
the transitional probabilities of Equation |8] are relevant for POAM analysis. I derive these equations and the variables contained 
therein in Appendix ICJ 

The third illustrative form of Iz is 



Iz ^ Y '"^»'.™ ^ ^"^ ['^^^ \ '^^^ n X R p(Q)/(R) = Yj "^Pm,m + 27r fd^nnpiO) X fd-RRp"{R) 



(10) 



where /?"(R) is the probability that a photon can pass through the aperture within d^R of R, and p(Q) is the probability that a photon 

arose from the celestial sphere within d^Q. of fi. The source-dependent term is expressed as Q x R operating on the probabilities or 
the cross product of the dipole moments of the probabilities. I derive these equations and the variables contained therein in Appendix 

o 

For an on-axis point source (p = 0), the POAM and torque spectra in the image plane are identical to the source-independent 
POAM spectra within the aperture, or ^,„,„, - r„,_,„ - p"„,„. Even if the object under observation is not an on-axis point source, the 
ensemble of ;?JJ, „, and the POAM quantities formed from them represent reasonable source-independent metrics. 

The source-dependent terms, on the other hand, are dipole moments with + 1 selection rules, which means that they are iden- 
tically zero on-axis and their effects are relatively small off axis. Elias (2008 ) called these terms "pointing" POAM or "structure" 
POAM. The source structure cannot be easily disentangled from the effects of propagation media and instrumentation. They are 
also zero when there is only a single non-zero Sk{R). 
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5. Interferometers 

lEliasI (I2OO8I) derived POAM correlations and rancors for a single-baseline optical interferometer They depend on baseline length, 
delay, telescope aberrations, etc. Since he was considering only a single observation with a single pair of telescopes and integrating 
over the image plane, employing the baseline vector instead of the two telescope vectors is acceptable. For this simple situation, the 
total POAM and torque can be set to zero because the synthetic aperture origin can be always placed along the line containing the 
single correlation measurement of two telescopes. 

In this section, I derive image-plane POAM quantities for an interferometer. There are slight differences between radio and op- 
tical interferometry. In the radio case, the electric fields between pairs of telescopes are multiplied and averaged. In the optical case, 
electric fields are summed, squared, and averaged. Mathematically, the same zero-spacing fluxes and visibilities can be obtained in 
both the radio and optical. I assume that the fringes are tracked well enough to avoid any delay dependence. 

Consider a single telescope behind an opaque aperture that contains a finite number of imperfect pinholes with amplitude and 
phase errors. The resulting "dirty" image is the convolution of the perfect diffraction-limited image and the Fourier transform of the 
imperfect pinhole pattern, which is equivalent to the "optomechanical" Fourier transform of the unnormalized visibilities from all 
pinhole pairs. An interferometer works in a similar manner Multiple observations with an array of small telescopes, representing 
the imperfect pinholes, form the "synthesized" aperture of a single large telescope (cf. Figure |2]l. Unnormalized visibilities of all 
telescope pairs are then mathematically Fourier transformed to create the dirty image. The sky-dependent response of the individual 
telescopes defines the image FOV and their imperfections also affect the images (cf. Sections [5. 2l andl6Tl. 

5.1. Initial Mathematics 

The single-telescope POAM mathematics of Section 143] can be used directly with interferometers, taking into account the discrete 
sampling of the aperture. For the sake of illustration, however, I rewrite them in terms of sampled unnormalized visibilities, which 
are the standard interferometer observables. The interferometer intensity response is 






E(n ■ t) 



^ Sr' e'^"^-^ .9(R')£(R;f) \ 
= { cfR' ^ cfRe-J^"'^^-'^^^ i(R)/(R)T(R,R) = { cFr' frf^/Je"^^''*"^'" ^(R ,R) , (11a) 

where £(R ; t) is the electric field in the synthesized aperture, i(R ) is the telescope-based gain function of the synthesized aperture 
(analogous to the aperture function of a single telescope, cf. Section©, ;F(R , R) is the uncalibrated unnormalized visibility, and 

r(R,R) = m'-R) = /^£(R;0£'(R;0) = (^ ^ cf-Q." e^-"^'^ E(il'-t) f d^ne-J^''^-^E''(Q,;t)\ 



is the true unnormalized visibihty under SAA. Note that the image-plane intensity of Equation II lal is expressed in terms of two 

aperture-plane integrals, as opposed to the standard single integral over baseline b = R - R (the separation of two telescopes), 
because POAM quantities depend on telescope position vectors not baseline vectors. In Appendix |D] I prove that this equation 
can be converted to the baseline form used for non-POAM analysis. I assume that a single moving baseline produces all of the 
unnormalized visibilities. All formulae, however, can easily be generalized to multiple moving baselines. 

Expanding the exponential functions in Equation ! llal in terms of Bessel functions, the (m,m)''' POAM autocorrelation density 
becomes 

B„,,„,(p') = fd'-R' fd^RJ„,(2^R'p')J„,(27:Rp')e'^"'^'^'"^^f(R,R). (12a) 

Integrating over the image plane, I obtain the (m, m)''' POAM autocorrelation 

B„,,„, = lim 27r dp' p' B,n,n{p') ^ 2n \ dRR — dt/^' — dif, e-J""-'"'-'^^ r(r,R) , (12b) 

Prav^<» Jo Jo 27r Jo 2.7T Jo 

where r = (R cos t//',R sin iff') and /?,„, is the radius of the synthesized aperture, both in units of wavelength. These equations are de- 
rived in AppendixlDJ As in the single-telescope case (cf. Section l4TTT i. Equation ! 12bl can be used to create all of the POAM quantities 
defined in Section[3]as well as the illustrative forms of Section l4r2l (cf. AppendixiDb. There are, however, two complications. 

The only uncalibrated unnormalized visibilities that contribute to the m''' POAM state autocorrelations are those which come 
from pairs of telescopes in the same aperture ring R, i.e., pairs of telescopes that are the same distance from the aperture reference 

point (/? = |r| = |R I = |R|, cf. Figure |2]l. This result is not surprising given that POAM quantities are defined in rings. As a matter 
of fact, it is possible to rewrite Equation ! 12b! as azimuthal Fourier series components of the azimuthal convolution of uncalibrated 
aperture electric fields integrated over radius. Unfortunately, real interferometers do not have telescopes arranged in this manner, 
which means that the true unnormalized visibilities and telescope-based gains must somehow be interpolated onto a polar grid. 



N.M. Elias II: Photon orbital angular momentum and torque metrics for single telescopes and interferometers 

Single telescopes obtain data using an entire aperture. Source-independent POAM quantities calculated from an on-axis point 
source calibrator observation can be used to judge the quality of separate science target observations if the atmospheric statistics are 
X! consistent and the source structure does not extend too far from the FOV center This strategy does not work for interferometers. 
They obtain data through a sampled synthetic aperture, so the sample coverage for an on-axis point source calibrator and a science 
target will be different. 

Determining the optimum strategy to overcome these complications requires a significant amount of effort. Such work is beyond 
the scope of this paper, but here I present two possible candidates that act as starting points for future research. 

When the ungridded discrete Fourier transform (DFT) of science target uncalibrated unnormalized visibilities is calculated (no 
additional processing; e.g., CLEAN), the resulting image is corrupted by incomplete sampling of the synthetic aperture and gain 
errors. If the inverse DFT (IDFT) of the corrupted image is calculated on a polar grid, it effectively interpolates the uncalibrated 
unnormalized visibilities so that they can be used directly in Equation ll2al to form the POAM quantities of Sections[3]and|4] 

Many interferometry imaging-processing techniques iterativel y solve for sampled true unnormalized visibilities and telescope- 
based gains while improving the image model (Ra u et al. ..2009: Raul 1201 Oh . If the synthetic aperture is sampled densely enough, 
the true unnormalized visibilities and telescope-based gains can be interpolated onto a uniform polar grid so that the desired POAM 
quantities of Sections[3]and|4]can be determined. In AppendixlE] I show that the interpolation kernel must be azimuthally symmetric. 



or TCI 



•R-R 



), in order not to further modulate the POAM spectrum. 

These two strategies yield different results. The DFT/IDFT interpolation method includes the effects of both incomplete sam- 
pling and gain errors. The azimuthally symmetric interpolation kernel method, on the other hand, includes only the effects of gain 
errors if the processing successfully removes artifacts due to imperfect synthesized aperture sampling. The DFT/IDFT interpola- 
tion method includes pointing/structure POAM which cannot be easily disentangled from source-independent terms (if the object 
under observation is an on-axis point source, there is no pointing/structure POAM). Conversely, the azimuthally symmetric in- 
terpolation kernel method estimates the aperture functions, which means that pointing/structure POAM can be disentangled from 
source-independent terms. 

5.2. Telescope Apodization 



The true unnormalized visibility (Equation [TTb} for a single baseline, modified by the apodization of the individual telescopes, is 

:r(R,R) = r£/2Qe>2<5-5).H j?i'(n)j?[*(n) B{n) = f d^ne^^'"^^-^>-^'P(n)B{n), (U) 

where J?I'(n) and J?l(Q) are the complex electric-field sky-dependent gains of two telescopes at points R and R in the observation 

plane, and !P(i2) is their combined power sky-dependent gain. This equation can easily be generalized for multiple moving baselines. 
If the electric-field gains are identical, the power gain is real. If the converse is true, the power gain is complex. A complex power 

gain represents a non-Hermitian calibration error, ;F*(R , R) ^ !F(R, R ), which leads to a complex output image (cf. EquationfT4l). 
Substituting EquationfTSJinto Equations[TTh-b. I obtain a modified version of Equation |4a] 

B(n) = fd^np(nM)P(n)B{h.). (i4) 

When expanded into POAM components, this EquationfT4lbecomes 

Bin) ^ ^ J] Bp,,ip')eJ^P-''^'^' , (15a) 



where the POAM correlations 

CO oo 

Bn 



oo CO 



m=—oo n=-oo *^^' lk=-oo l=-oo 



B,„,„ip) (15b) 



have an extra function 

n-/-,„+„(p) = ^ ( d4>e-J'-''-'-"'^"^'^PiQ.). (15c) 



The interferometric PSF input/output (separate) gain is 

P;,'f'(p'^p) = D;,\p',p) D-;*(p',p) , (16a) 



where 



^«>''P) = ^r d4.e^"'^D„,(p',n) (16b) 
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is the interferometric DF input/output gain. These gains are defined in Tables 2 and 4 of lEliasI (l2008l) . I derive Equations [T5b-c in 
Appendix|F] 

From Equation I15bl I see that sky-dependent gains do indeed modulate POAM. To understand these effects more clearly, I 
choose a simple case where the interferometric synthetic aperture is fully sampled with no amplitude or phase eiTors, which means 
that 

Ptt^P'^P) ^ P],T'^'^P)^Kph, (17a) 

(lEliasll2QQ8b . Equation ll5bl then becomes 



CO CO 



Bp,,(p') ^ Yj Yj^"" \ dpp[Pp'^-'(p',p)'Pp-,-„,^„(p)]B,„Jp). (17b) 

ttt — — rv> n — - — rv^ *-' V7 



m=-co «=-oo 



Now only the sky-dependent gains distribute the input POAM correlation densities to multiple output POAM correlation densities. 
The index of the extra function consists of the difference of image plane and celestial sphere index differences. To further verify that 

these mathematics are correct, I let the sky-dependent gain exhibit only radial apodization, or 'P(Q) — > 'P(p)- The new extra function 

"Pp-q-m+nip) ^ foip) Sg,p-,„+„ (18a) 

produces a one-to-one correspondence between input and output POAM correlation densities 



B„,,„(p') ^ 2n \ dpp y PJ,j;Z!!XT"(p''P) 
Jo ,r^^ 



p,p-m+n 
lP=- 



■Po(p)B,n,„(p), (18b) 



as expected. 

The sky-dependent gain for each antenna can be determined via holography, i.e., raster scans around a blight point source 
normalized (amplitude and phase) to a reference antenna (cf. Section |6]l. The Fourier transforms of the sky-dependent gains and 
their inverses are 



a'(r') = fd^Qe-^^"-^Jl'(n) £ ^'(Q) = f d^ r' e'^" '^ a' (r' 



) 

(19) 



where a'(r ) and a(r) are the complex holography functions representing aperture imperfections projected back to planes in front 

of the telescopes. The aperture coordinates are r = (r'cos;^'',/ sin;^'') and r = (rcos;t',rsin;t^). These functions are conceptually 
identical to the single-telescope aperture function (cf. Section 14.11 1 and interferometer synthesized aperture function (cf. Section 
15.1b . The POAM components of the inverse transforms are 

dr' r' Jk{2nr'p)a[{r') and J[t(p) ^ f 2n \ drr Jk{2nrp)ak{r) , (20a) 

Jo 

where 

-2n 






k=-oo 

^kip) = — dcl)e-i'"^J[{n) « miO) = y ^k{p)e^''^ (20b) 
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2^ Jo M^ 



akir) ^ — \ dxe-J'^air) « a(r) = 2j «*W^ • (20c) 

k=-co 

->' -4 

Note that when the apertures are azimuthally symmetric, or fl'(r ) — » a'(r') - a'^ir') and fl(r) — » a(r) - ao(r), only the ^'q(p) and 
^o(p) terms are non-zero. Further, if both antennas of a baseline are azimuthally symmetric, it follows that their power pattern is 

also azimuthally symmetric P(Q.) —> "Pip) = 'Poip) = ^o(p) ■^o(p) and does not redistribute POAM states. 
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6. EVLA Holography and POAM 

Recent K-band (^ 24 GHz) holography observations of EVLA antennas were processed during commissioning^ (lBrentiensll2011h . 
The target is a bright point source. Three of the antennas tracked the target and were used as amplitude and phase references. The 
rest of the antennas - i.e., those under test - executed raster scans. The resulting data were flagged, calibrated, and averaged before 
calculating the Fourier transforms on Cartesian and polar output coordinates. Calibration included removing the effects of pointing, 
focus, and subreflector rotation offset, so the holography should ~ represent the imperfections of the antenna surfaces. 

The amplitude and phase responses (Cartesian coordinates) are shown in Figures |3] and |4] When calculating the Fourier trans- 
forms, I oversampled the output by a factor of six for smoother interpolation. The circular shape of the dishes is clearly visible. 
The amplitude responses exhibit an opaque region near the center and four orthogonal "spokes" produced by the subreflector and 
its supports. There are also clear indications of both large and small spatial scale reflectivity features. The phase responses show 
random phases in the opaque regions, as expected. The locations of the small scale phase features (mottled patterns) are consistent 
with individual misaligned panels on the reflector surface. They tend to appear in groups. 

In Figure|5] I display the sky-independent image-plane POAM probability spectra (cotTesponding to the aperture POAM proba- 
biUty spectra in Equations|8]|9al andfTOli as well as the total torques. There are no sky-dependent quantities (i.e., no pointing/structure 
POAM) because the object under observation is an on-axis point source. To determine the sky-independent image-plane probabil- 
ities for each telescope, I calculate the azimuthal Fourier series components of its aperture function versus radius flt(r), form the 

squared magnitude of each Fourier component versus radius \aic{r)\-, sum each squared magnitude over radius to obtain la^l^ (times 
27r r Ar, where Ar is the radial size of the aperture element), and normalize the la/tl^ to the sum of the \ak\^ over k. The range of 
POAM components is limited to + 15, which is w the Nyquist limit at the largest radius (I do not calculate higher POAM components 
even though they are available because of oversampling). 

When a perfect telescope observes an on-axis point source only the n = component of the image-plane POAM spectrum is 
non-zero. Conversely, when an imperfect telescope observes an on-axis point source the n -Q component is reduced and the other 
components become non-zero. Most of the n + Q probabilities are 1% or less. A few of the n — +\ probabilities are almost an 
order of magnitude larger, which could be caused by feed position errors or uncalibrated pointing errors. Summing over the n + Q 
probabilities, we find that ~ 10% of all photons are "torqued" away from the « = state. Also, note that most telescopes exhibit a 
non-zero total torque but some (e.g., telescopes 24 and 28) exhibit w zero total torque (i.e., oi symmetric POAM spectra). 

7. Conclusion 

With previously defined concepts and calculi (' Eliasll2008l) . I presented generic expressions for POAM spectra, total POAM, torque 
spectra, and total torque in the image plane. I extended these functional forms to describe the specific POAM behavior of both 
single telescopes and interferometers. These POAM quantities make excellent metrics, complimenting Zemike polynomials, for 
describing the response of astronomical insttTiments. Real holography measurements of EVLA antennas demonstrated their utility. 
Now that POAM metrics have been derived, it is incumbent on the author to make them available to the astronomical commu- 
nity in an imaging package. In the future, I plan to extend them to handle spin-polarized and non-flat spectrum sources. Possible 
future studies include forming POAM quantities using real interferometry visibility data, modeling instrumental aberrations and 
turbulence of the tropos phere/ionosphere in ter ms of POAM, POAM-based imaging algorithms and constraints, POAM-based su- 
per resolution imaging (iTamburini et al. I [20061) without interpolati on or extrapolation, and POAM observations of astrophysically 
important sources such as masers and black holes (JHarwit 2003; T amburini et al. 11201 lb . 

Acknowledgements. NME2 would like to thank Dr. Sanjay Bhatnagar for fruitful discussions and advice, and Drs. Michiel Brentjens, Richard A. Perley, and Bryan 
J. Butler for providing cahbrated EVLA holography data. 
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Fig. 1. Schematic diagram of a single telescope looking at a source on the celestial sphere. The coordinates are indicated and 
described in the text of Section|4] The electric fields (celestial sphere and image plane) and aperture function are also presented. 
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Fig. 2. Schematic diagram of an interferometer looking at a source on the celestial sphere. The coordinates are indicated and 
described in the text of Sections |4] and |5] The electric fields (celestial sphere and image plane) and aperture function are also 
presented. The aperture function refers to the synthesized aperture, not to the individual telescope (or pinhole) apertures that sample 
it. 
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Fig. 3. Holographic amplitude measurements of the EVLA antennas. The blue and red colors represent the low and high reflectivity 
regions. The coordinate units are wavelengths. 
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Fig. 4. Holographic phase measurements of the EVLA antennas. The yellow colors represent raised regions with respect to the 
fiducial dish surfaces. The RMS panel deviations are s; 200 yum. The coordinate units are wavelengths. 
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Fig. 5. Holographic probability spectra of the EVLA antennas. The ordinates are log (base 10) probability/100% and the abscissae 
are POAM states. The total torque is written on each subplot. 
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Appendix A: POAM from the Celestial Sphere 

lEliasI (I2OO8I) tacitly assumed that the total POAM on the celestial sphere was zero for SAA sources and non-zero for non-SAA 
sources. In this appendix, I formally prove those statements for a single telescope, but the derivations apply to all optical systems, 
including interferometers. 

I consider the SAA case first. If the instrument does not modulate the POAM spectrum, the aperture function in Equation l5fl 

becomes ^(R) = so(R). The integral function in Equationl5e]then simplifies to .T,, Jo', o) = J^p,p(p',p)Sq^p. When this result is substi- 
tuted into Equation|5c] the (m, m)''' POAM sensitivity is no longer a function of (f>, or Pm.mip', il) -^ Pm,m(p',p) - \j^m,m(p' , p)\ ■ Upon 
inspection of Equation l5el I find that J'-m.-mip' , p) - J^m.m(p',p), which means that f _,„,_„,(p',p) = Pm,m(p' iP)- Since P,„,m(p',p) 
is even in m, the total image-plane POAM of Equation [7] is identically zero. Because I initially assumed that the instrument does 
not modulate POAM, it follows that the total POAM on the celestial sphere must be zero as well. Further, the total POAM on the 
celestial sphere for SAA sources must always be zero because the behavior of the source must be completely independent of the 
behavior of the instrument. Q.E.D. 

I now consider the non-SAA case using a completely spatially correlated source. The spatial and temporal parts of the electric 

field on the celestial sphere factor into separate functions, E{Q,; f) - E{Q,) f{t), where f{t) is a random complex function. The 
electric field in the image plane (Equation |4a|i becomes 



E(n ; t) 



I 



cf-Q.D{Q. ,n)E(n;t) 



fit) =£(Q)/(0. 



fd^Q.D(n,Q:)E(n) 

Expanding the DF into sensitivities (Equationl5db yields 

00 00 CO r ^ 

£(n;0 = 2 £„,(p';Oe>'"*' = Ji [E,n(p') fit)] eJ""^' = Xi / ^'s^^m(p',i2)£(n)/(f) 
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Jmf 



(A.la) 



(A. lb) 



Apart from the separation of the spatial and temporal components of the electric field and the POAM states, these functional forms 
are identical to the SAA ones. 

With these POAM states, I form the image-plane POAM autocorrelation densities (Equation|5a]i and express the PSF sensitivities 
in terms of the integral functions (Equation l5ct 



B,„,Jp') = i^\E,„(p';tf 



fd^QD,„(p', 



£2) E(n) 
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(fQ.J,„,k{p',p)e-'''*E{m 



(A.2a) 



where F - < \\f{ty^ > is the RMS of the random complex function. Again I assume that the instrument does not modulate POAM, 
or J^m,k{p\p) - Sfm,m{p' ,p) 5k,nu thereby collapsing the sum in this equation. Combining this fact with integration over radius p' in 
the image plane, I obtain the total (m, m)'* POAM autocorrelation in the image plane 
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This equation shows that B-,„.-,„ ^ B,„.,„ because E-„,{p) + E^ip) in general, which means that the total image-plane POAM can 
be non-zero (Equation [3]l. The total image-plane POAM is the same as the total source POAM because the instrument does not 
modulate POAM, so if the former is non-zero the latter must be non-zero as well. Q.E.D. 

In the previous two proofs, I solved for the POAM autocorrelations in the image plane and inferred the form for the POAM 
autocorrelations on the celestial sphere. Here I provide a bonus proof employing some mathematics that are not found elsewhere in 
this paper. I only use the POAM autocorrelations on the celestial sphere and I express them in terms of azimuthal correlations of 
electric fields. This proof can be used to describe the POAM behavior of SAA and non-SAA sources. 

The (m, m)''' POAM autocorrelation on the celestial sphere is 
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Jo 
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is the temporal correlation of the electric fields on the celestial sphere between azimuths 4>' and (p integrated over radius p, and 

B(p,0',0) = r-E{p,4>'-t)E*{p,(l,-ty) = /^[£.(p,<^';f) + ;£,(p,0';f)][£.(p,0;O-;£,(p,0;O] 



^-E,(p, 0'; t) E,(p, 0; f)\ + l^Eiip, ip'; t) Ei(p, 0; t) 



+ J 



^E,(p, iP'; t) Eiip, cf,; t)\ + l^Eiip, cf>'; t) EAp, <p\ t) 



(A.3c) 



= [Brrip, 4>\ <p) + Biiip, <f>', <f>)] + j [-B,,(p, 4>', 4>) + BiAp, <p', 0)] = Br(p, <p' , <p) + ;B,(p, <p\ <p) 

is the temporal correlation of the electric fields on the celestial sphere between azimuths (p' and at radius p. In general, 

and B{p, (p' , <p) are complex numbers written in terms of complex electric fields. Similarly, the (-m, -m)''' POAM autocoiTelation on 

the celestial sphere is 
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j^ dppl^ |£_,„(p; t)A ^ l-j^Jd<p' e-J""^' ^^J dcp e>'^ B*(^', 0) , (A.3d) 



(A.4a) 



B-m,-m = 27r I dppB-„,_^,„(p) = 2n I dpp l-\E-,„(p;t)\ 

It is obtained by replacing m — » -m and exchanging 0' <-^ in Equation lA.3al 
As stated elsewhere, fe = when B-„,.-,„ - B„,„ for all m. I define the quantity 

AB,„ = B,„,,„ - B_,„,_„, = ;2:1 f rff g-^'^'^i f dc(>e^""^ Bi(cf,',<f,), 
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where 
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B!(cf,',<p) = 27r I dppBi(p,4>',4>) = -27r ( dp pB,,{p,(t>\<p) + 27r ( dpp BUp,<f>',ip) = -B,-i(<f>',<f>) + B i,(<f>' , cf>) . (A.4b) 
Jo Jo Jo 

For an SAA source Equations rA.3cl and [Ob] lead to B ri((f>' , (p) - Bir((f>',(f>), Bi((p',4>) - 0, AB^ - for all m, and thus Iz - 0. In 
other words, if the complex electric fields are spatially uncorrected, their real and imaginary parts are spatially uncorrected as well. 
Conversely, for a non-SAA source I find that Bri((f>', 4>) + Bir{4>', (p), Bi((p', cp) + 0, KB^ + for all m, and thus Iz + 0. Q.E.D. 

Appendix B: Telescope POAM Autocorrelations 

The {m,m)'^ single-telescope PSF sensitivity (Equation [Sc]) is the squared magnitude of the m'^ single-telescope DF sensitivity, or 
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The m single-telescope DF sensitivity is simply the azimuthal Fourier component in the image plane 
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The quantity in square brackets is j " Jm(2nRp') e ^'"'^ . When the aperture function is expanded into azimuthal Fourier components, 
I obtain 
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The quantity in square brackets is f^'^ J p-„{2nRp) g-'^''"'"''^. After replacing p with m - k and rearranging, the result is 
D,„(p',£l) = j-'" 2 / 27r dRRJJ2nRp')s„-k(R)Jk(2nRp) g-'*^ = f"" ^ f J,„,k(p' ,p)e-^'"^ . 
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The squared magnitude of this equation is identical to Equation|6b] which can be substituted into Equation |5a] to give the complete 
(m, m)'^ single-telescope POAM autocorrelation densities. Q.E.D. 



Appendix C: Illustrative Forms 

To derive the first illustrative form, I change the indices k —> m - p and I —> m - q of Equations |6b] and [T] The resulting total POAM 
kernel becomes 
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When Equations lC.lb -b are substituted back into Equation|7] I obtain 
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is the integrated and squared (m, m)''' component of the aperture function, the 
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are the aperture dipole-moment functions, and the 
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are the first-order source rancors (lEhasI EOOSl Equations 12a-b). Note that p was replaced with m so that Equations IC.2fa -c are 
consistent with the equations in Section l4~2l The total intensity in the image plane can be rewritten as 
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The quantity in the square brackets approaches 5(R - R ), which means that 
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where ^(R) - \s{R)\^. When this equation is substituted into Equation IC.2al I can define the aperture probabilities /?Jj, ,„ - S„,.m 
I S , the aperture transitional probability densities 75° ^j^j(/?) = «„,(/?) .s*j^j(/?) / S , and the celestial sphere transitional probability 
densities p„j,±i(p) - B„.„±i(p) / B. Keeping in mind that S_i(p) = S* j(p) and D-(p) - D\{p), the result can be rearranged to obtain 
Equation!] Q.E.D. 

To derive the second illustrative form, I rewrite the sum over transitional probabilities in the aperture in Equation |8]as 
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where R - {Rcosil/,Rsmt//) and r - (Rcost//' ,Rsmt//') (both aperture points are located in the same ring /? = |R| = |r|). The 
quantity in parentheses is 2;7r 6(il/' - t//), so 
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which are the normalized first-order aperture rancor gains. Substituting Equation IC .4b I into Equation IC.2cl Equation IC.2al can be 
rearranged to form Equation |9a] Q.E.D. 

To derive the third illustrative form, the imaginary part of the integrand of Eauationl9alcan be simphfied to 
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When this expression is substituted back into Equation|9al Equation[TO]is obtained. Q.E.D. 



(C.5) 
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Appendix D: Interferometer POAM Autocorrelations 

To prove that Equation lllal can be expressed in terms of a single integral over baseline for standard non-POAM interferometry 

— > — >' — > 
analysis, I substitute b = R - R to obtain 

B(Q)= fd^R fd^be-J^''^-^s(R+h)s'-(R)rih)^ f d^be'^^"^-^ f d^Rs(R + h)s\R) Tih) = f d^be-^^"^'^ SihWih), (D.l) 

where 5 (b) is the baseline-based gain function. The +R shift in the b integral can be ignored. The integral inside the square brackets 

must be the baseline-based gain function. Each point in the aperture function s(R) corresponds to the position of a telescope. The 
aperture function appears twice in the integral because there are two telescopes in a baseline. The two telescopes are separated by 
the baseline. Q.E.D. 

To derive the (m, m)''' interferometer POAM autocorrelation density, I substitute 



into Equation ! 11 al and find that 
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Each m = n term is identical to Equation 1 12al When Equation ! 12al is integrated over radius in the image plane, this integral 
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appears as one of the factors in the result. This Dirac delta function collapses the R' integral, leaving Equation ! 12bl Q.E.D. 

To prove that the interferometer POAM autocorrelations of Equation ! 1 2b! can be converted the single-telescope illustrative forms 
(cf. Section l4r2b . I form the total POAM in the image plane 
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h ^ 2_j '^P'"'"' ^ "fi Zj '"^"'."i ^ "fi^^ I dRR— I dil/'— I di// 2_j 

m--co ;;?=-oo ^" ^" ^" Lm=-oo 

The quantity in the square brackets is 

CO 

^ me-'"'^''''-'^'^ = j In 6' (iff' -if/), 

m=-co 

which is proportional to the derivative of another delta function identity. Since 

I dxg(x)6'(x-xo) = -8'(x)\^^^^ , 



-jm(ili'-ip) 



T{r,R). 



(D.3a) 



(D.3b) 



(D.3c) 



Equation !D.3a! then becomes 

I2 = -j fd^R '^^'"'^^ /B = -;• fd^R -4-^ s*(R)/S + f d^Q f d^R [2nRpsm((p- dr)] Bin)S(R)/B 

J dip' J dt//' J J 



<P'=<P 



£ m<„, +27rJ^2^J 



d^R 



QXR 



'/''='A 
p(Sl)p\R), 



(D.3d) 



where the required variables and functions are defined in Appendix ICl This equation is identical to the third illustrative form. The 
third illustrative form was derived from the other two in Appendix|C] so there is no need to rederive them here. Q.E.D. 
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Appendix E: Interpolation with No POAM Modulation 

In this appendix, I assume aperture function interpolation of the form 
i(R) ^ s(R) = f d^'R'KiR,^) sm , 



(E.l) 



where '7C(R, "R) is the interpolation kernel, R = (R cos i//, R sin tp) is the interpolated aperture coordinate, and R-CR cos *?, R, sin *P) 
is the real aperture coordinate. Both coordinates are in units of wavelength. This equation can be expanded into POAM components 



S(R) = y UKJe''"^ o UR) = — di/,e-J"'^S(R) = V 27r dRR%;"\R,R) s,„{R) , 

„t-L 27rJo „^^ Jo 

where Z?/,,,,,, is the radial integration limit (Ri^i for single telescopes or /?,,„ for interferometers), 

'K(R,'R) = y y 'K-"'iR,'R)e--'"'^e-'"'^ o 'K;"'{R,'R) = — di/ie-''"^ — 

„ti„,t-L 27rJo 27rJo 

is the double Fourier series expansion of the interpolation kernel, and 



-.infl. f d'V e^"'^ 'K(R,R) 



(E.2a) 



(E.2b) 



(E.2c) 



is the POAM expansion of the non-interpolated aperture function. All of these equations have been derived using POAM calculi 
(lEUasll2008l) . 

By definition, a system that doesn't modulate the POAM spectrum means that a single input POAM state gives rise only to the 
same output POAM state ( Elias 2008 ), so the interpolation kernel POAM gain must be 'K^"'{R, R) - K„,{R, R) 6n,m- The interpolation 
kernels that produce such POAM gains consist of the family of functions that exhibit circular symmetry, or 



'K{R,R) = 'K( 



R-R 



) = "Ki ^R'^ +R^ -2RRco&{tf/ -"¥)) . 



To prove this assertion, I substitute this kernel into Equation IE. 2bl and obtain 

'K;'"(R,'R) = — r dxe-J'"^'K{jR^T^'^^2m^^I(ip^^)) — f 
[27r Jo ^ J [27r Jo 

which means that 

S„(R) = 27r d'R'RK„(R,R)s„(R). 



d)p 



,-j{n-m)il/ 



K,„(R,R)6„,, 



(E.3) 



(E.4a) 



(E.4b) 



I employed a change of integration variable ^ = i// - x and eliminated i// from both x integration limits. Q.E.D. 

This type of circularly symmetric interpolation kernel is often used in interferometiic imaging, which is indeed fortunate. 
According to Equation IE. 4bi it appears that the s„{R) must be calculated in order to determine the Sn(R). This procedure is unnec- 
essarily complicated. If one substitutes Eauation lE.2cl into Equation IE. 4b I one finds that 



S„(R) = fd^K [K„iR,'R)e-'"''] s(K) . 



(E.5) 



In other words, the n''' interpolated POAM gain can be calculated directly from the continuous or discrete aperture function and the 
kernel contained in the square brackets. 

Appendix F: Telescope Apodization in Interferometers 

The interferometric PSF in Equation fl4lcan be expanded as a interferometric DF dual azimuthal Fourier series squared, or 



p(n , n) 



D(n , a) 



OO CO 



CO CO oo oo 



p=—o3 k=-oo p=-co q=-oo k = -oo /=-co 

where D~"(p\p) is defined in Equation ll6bl Similarly, the source intensity can be expanded into its POAM correlation expansion 



CO CO 



Bin) = 2 Z 5«v.(p)e*'-">^ 



(F.2) 



m——co ?7=— oo 



When these equations are substituted into Equation [14] the indices rearranged, and the integral over (p performed. Equations [T5b-c 
result. Q.E.D. 
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